Analysis of "Gauge Modes" in Linearized 
Relativity 



Richard A. Matzner, Mijan Huq, Alonso Botero, Dae II Choi, Ullar 
Kask, Juan Lara, Steven Liebling, David Neilsen, Premana Premadi 
and Deirdre Shoemaker. 

Department of Physics and Center for Relativity, 
The University of Texas at Austin. Tx 78712-1081 

Abstract. By writing the complete set of 3 + 1 (ADM) equations for linearized 
waves, we are able to demonstrate the properties of the initial data and of the 
evolution of a wave problem set by Alcubierre and Schutz. We show that the 
gauge modes and constraint error modes arise in a straightforward way in the 
analysis, and are of a form which will be controlled in any well specified convergent 
computational discretization of the differential equations. 
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1. Introduction 



The 3 + 1 (ADM form of the Einstein equations for vacuum may be written: 



R + K 2 - K cd K' 



cd 







(1) 

(2) 
(3) 
(4) 



V a K a b -V b K = 
d 9ab = ~2aK ab 

d Q K ab = -V a V b a + a[R ab - 2K al K\ + KK ab ] 



where the derivative [2] d Q is defined: 



d — do + £p- 



(5) 



Here a is the lapse function, which relates the coordinate time between time 
instants to the proper time interval, (3 l is a shift vector describing changes in the 
coordinization at one time and £p is the Lie derivative along (3 % . Because they 
are kinematical choices, we can, for some (perhaps small) interval, set a and (3 l 



posed the problem we now define. We use the overdot and d a interchangeably below 
to indicate the time derivative. 

In Eqs. (1-4), all of the geometrical variables are 3-tensors. g ab is the spatial 
3 metric; K ab is a spatial tensor (with a = 1, (3 % = 0, equal to — \g a b)', K is the 
3-trace of K ab ; R ab is the 3-Ricci tensor obtained from g ab , and R is the trace of 
R ab . Linearizing the equation around flat space (g ab = 5 ab + h ab \ K ab "small") and 

o o 

introducing the traceless variables K a b= K ab — -^5 ab K, and h a b= h ab — ^h^d 5 ab , we 
have 



arbitrarily. Here we take a — 1, /?* = 0, following Alcubierre and Schutz, || who 



R = 

2 

K a b) a = ~K, b 



(6) 
(7) 
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h = -2K 



(9) 
(10) 

(11) 
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mn,ab 



(12) 



where the comma now denotes the partial derivative. 

Formally, Einstein equations are a constrained hyperbolic system. Solution of 
the constraint equations (1 - 2) [or their linearized version (6 - 7)] will be preserved 
analytically in the evolution of the dynamical equations (3 - 4) [or their linearization] . 
By setting a = = 0, we choose not to consider the admissible possibility of 
perturbations in a and f3 l of the same order as h. This is a coordinate choice (gauge 
choice) . Further, there is the possibility to specify the gauge more completely at the 
perturbation level in a way that does not change a, (3 l . Alcubierre and Schutz [3], for 
instance, choose a transverse traceless gauge, S mn h mn = and 5 mn d m h n i = 0. Such 
a choice can be accomplished by a coordinate transformation at the original time; 
analytically this choice is preserved with a = 1, (3 l = in the evolution. Gauges like 
this represent additional restrictions on the behavior of the solution; however they 
are not necessary to have a physically meaningful evolution. We shall see below that 
even in this linearized case, setting such gauges in computation is delicate. 

2. "Gauge Modes" 

Alcubierre and Schutz wrote the linearized evolution equations (8 - 11) in second 
order form without separating out the trace: 



With this form, Alcubierre and Schutz discovered secular, non-propagating, 
modes in their numerical implementation. In their attempt to evolve transverse 
traceless waves they find offending extra modes that violate the TT-gauge condition. 
Some of these modes appear to violate the constraints. By investigation of the 
properties of Eq. (13), Alcubierre and Schutz were able to show analytically the 
existence of modes similar to those seen numerically. 



{-dl + V 2 )hij = 5 mn (d n dih jm + d n djh im - didjh mm ). 



(13) 
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3. Preliminaries 



From Eq. (12) for the linearized Ricci tensor, it is easy to obtain the linearized 
3-scalar R: 

R=-S/ 2 h + h mn , mn (14) 

and the linearized Hamiltonian constraint (6) sets this to zero. 
A basic decomposition which we employ is: 

h - h TT 4- h L -+- -r) h* (iM 
lb mn "'mn ' IL mn ~ g u rnn lb ■ \- L,J J 

Here ti^ n is transverse: 

h TTmn , n = 0, (16) 
and traceless: 

h TT \ = 0. (17) 
The quantity h^ n is longitudinal, and can be written: 

for some vector v m . Notice that h^ n is not tracefree: 

h L m m = 2v m > m (19) 

(The construction (18 - 19) for h L mn is similar to that given by York|J except that 
we do not subtract the trace from h L mn ). Additionally we posit an independent 

o 

contribution to the trace, h*. At this stage we distinguish h mn , which need not be 
transverse, from h TT mn , which is transverse. 

Firstly, notice that h^ n annihilates the scalar curvature. This is not surprising 
because such longitudinal components arise from classic "gauge terms" as in 
Eq. (18). 

Secondly, notice that if we set (only) nonvanishing data, R = becomes 

V 2 h* = 0. (20) 

In this study we are concerned with computational relativity, for "physically 
realistic" disturbances. In particular, we will expect fall-off "at infinity" (this can be 

4 



accomplished by a mixed outer boundary condition on the computational domain), 
and we exclude singularities in the interior of the domain. The solution to Eq. (20) 
then is h — 0. Hence, we expect the trace, h, to arise only in the context of 
longitudinal components, h^ n , as in Eq. (19). 



4. Linearized Data Setting and Evolution 

We continue the Alcubierre and Schutz analysis of the linearized case, but take more 
careful note of the constraint equations in the analysis. We will also give a complete 
data-setting analysis. 

An immediate result from the Hamiltonian constraint R = and from the second 
order form of the linearized evolution equations (8 - 11) is: 

d 2 h = (!) (21) 

Hence the trace of the metric perturbation is solved by 

h = a{x l ) + 6(2;*)*, (22) 

if one enforces the Hamiltonian constraint to this order. Notice that since K = —\h, 
we have K = — \b(x % ). Notice also that because h can only arise from longitudinal 
modes, we have found Alcubierre and Schutz's growing modes: 

h L mn = A L mn {x l ) + tB L mn {x l ), (23) 



However, also note that B^ n is related to K mn and is restricted by the momentum 
constraint. 

o 

Similarly writing the second-order equation for hob, we find: 

(-d 2 + V 2 ) hij = 28 n d {t hi- \didjh (24) 

- l -(2d n d m °h mn - l -V 2 h)5 tr 

The last term (in parenthesis, proportional to |5y)in Eq. (24) is in fact part of the 
Ricci scalar, and again using the Hamiltonian constraint it equals —\V 2 h5if 

{-d 2 + V 2 ) hij= 2d n d (t hi- \ididj + 5^V 2 )h. (25) 

5 



Clearly, TT data annihilates the righthand side of Eq. (25), so h^ n satisfies the 
source free wave equation. It is also straightforward to verify using Eq. (18) that 
h I ^ nn annihilates the sum of the spatial derivatives in Eq. (25), so that 

d 2 h L mn = 0, (26) 

consistent with the behavior we have already found for the trace. Thus this system 
is simply hyperbolic in the notation of reference ||: disturbances travel at speed 
zero (h L mn ), or unity (h TT mn ). In order to proceed, we should also investigate the 
time derivative of the right side of (24). Let us proceed by taking note of one of 
the remarkable lessons learned from the hyperbolic analysis, for instance like that 
carried out in reference Q, that one should consider higher temporal derivatives of 
the evolution equations. For instance, an additional time derivative gives: 

{-d 2 + V 2 )K l3 =d l (K]) yn +d j (k?),n 

- -didjK - ^V 2 K. (27) 

By the momentum constraint, (7) the first two terms on the right are equal; 
each equal to \K^j. 
Thus: 

(-d 2 + V 2 ) Kij = di^K - -^KSij (28) 
= given function of spatial coordinates 
= -I(^.&-I^V 2 &). 

Hence; 

Kij = Kij (wave) + V~ 2 (didjK - ^V 2 K5 %J ). (29) 
Where V -2 is the Greens function with appropriate boundary conditions. 

.. o 

Furthermore, one can consider the temporal evolution of o Kij,i, the divergence 

o 

of Kij. 

By reorganizing Eq. (28) to have only the second time derivative on the left, 
and taking the divergence, we obtain 

<);, i< (30) 
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which in principle allows 

K ij tj = c%x j ) + d\x ] )t. (31) 

If we inspect Eq. (28) again, however, we find that solutions consist of solutions 
to the homogeneous wave equation (which, by Eq. (30) must be divergenceless) plus 
a "particular" solution which has no time dependence, because the right hand side 
is time independent. Hence di{x^) = in Eq. (30), and we may write 

Kij— K T J + V'^didjK - i(%V 2 K). (32) 

Consider an arbitrary longitudinal = w it j + Wj ti . A general 3-tensor of this 
form does not satisfy the momentum constraint j = SjKj: 



w, 



3 



+ w\j - 2wi' l yi = Wi'\i - w'\ij (33) 



Hence in order to satisfy the momentum constraint we must have 

"■i.,) J A l " 0. (34) 

The simplest solution is to take wi = ipj, the gradient of a scalar. Although 
more general solutions may be possible, but likely excluded by reasonable boundary 
conditions. With wi = ipj we have 

2V 2 ^ = K. (35) 



The corresponding K fj is 

2 1 
-SijV ip = 2ifj tij - - 



K % = - -(SijVV = - ~S tJ K. (36) 



5. Transverse Data 

As an aside, we state how one sets TT data. Obviously, if the data are set with 
a given (fixed) wave vector direction, one can algebraically make the tensors TT. 
Also, if one wishes to Fourier decompose the data, then each Fourier component can 
be made transverse. In general, one can use the following procedure. The idea is to 
pose an arbitrary traceless field C mn , and compute its divergence. 

Then, find the longitudinal component L mn = u> m ,n + w n ,m by solving 

C m % = V n (V m w n + V n w m ), (37) 
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i.e.: 

V n V n w m + V m (V • w) = C m \ n . (38) 

The operator on the left is easily derived from the minimization of 
V( m u>„)V (m iy n ' ) and so is strongly elliptic, guaranteeing the existence of unique 
solutions. If, for instance one postulates compact support for C mn , one can solve 
this equation using a scalar potential field <fi '■ w m = <p tm 

2V 2 0, m = C m \ n , (39) 

so that 

w m = <f>, m = V~ 2 {\c rn \ n ) ~ - + O(l). (40) 

The resulting longitudinal component L mn thus has behavior ~ 0(r -2 ), and the 
transverse data now is: C mn — L mn , and no longer has compact support. 

6. Complete Data Setting 

Set h^, h^ n = v m . n + Vn t m with appropriate locality (eg. compactness; at least 
"fall-off at infinity). Then the trace h is specified; cf. eqs. (19 - 20). 

Set K with appropriate locality. Then K^ n is set by Eq. (36) in terms of 
ip where 2V 2 ip = K and K^ m satisfies the momentum constraint. The Hamiltonian 
constraint is maintained by any K^ n , in particular one of form (36). 

7. Expected Numerical Behavior - Unconstrained Case 

o 

From Eq. (32) h ij satisfies 

h fj = h ij{to) ~ 2tV~ 2 {d i d j K - i%V 2 K) + /i5 T (wave). (41) 

Note that h ij (wave) is transverse because its time derivative —2K^ T is transverse. 

Since we have specified the intermediate variable ip by 2V 2 ip = K, we may also 
write: 



h ij —h 



TT 



(42) 



To proceed, let us suppose that we wish, as Alcubierre and Schutz did, to set purely 
TT data. 

In the Alcubierre-Schutz case, plane waves were set, so TT data can be 
algebraically enforced, and we are solving only 

Uh T Z = 0- (43) 

However, errors in setting TT data, can lead to a nonzero longitudinal part. Since 
we solve to zero the longitudinal component via an elliptic equation in a second 
order scheme we expect the longitudinal error to be second under small: 

l|A/iLll~ (^) HQ (44) 

where Ax is the discretization scale, and A is the typical scale of the TT part. Then, 
since longitudinal terms grow as t, one expects longitudinal contamination equal to 
the transverse signal at 

' AX )- 2 - (45) 



A v A 

For typical discretizations, (^f ) < 1CT 2 at worst, which suggests evolutions for 
times 

t ~ 10 4 A (46) 

before the longitudinal signals are comparable to the transverse signal, and another 
factor of order O before the longitudinal mode violates the linearization 

criterion ||/^ n || -C 1. Unfortunately, a poorly controlled computational scheme 
could put error (noise) in the longitudinal mode at short scales (so that ^ ~ 1). In 
that case one would expect to see the longitudinal mode grow to the amplitude of 
the transverse signal after only a few time steps; in a poorly designed differencing 
approach these secular zone-to zone oscillations would crash the program shortly 
thereafter. Notice that errors in the trace h will grow similarly to the h L mn modes 
and with potentially worse effect, because they violate the Hamiltonian constraint. 
Without a close inspection of the Alcubierre-Schutz code, we cannot comment on 
why they found such rapid growth of non- TT modes. 

In a recent paper 0, Gundlach and Pullin point out a mechanism for 
instabilities arising from the violation of the constraints in a free evolution. They 
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used perturbation analysis in double null coordinates on a Reissner-Nordstrom 
background, and found that a free evolution led to exponentially growing gauge 
violating modes. Their results can be taken to the flat space limit. In contrast to 
the Gundlach and Pullin result, our flat background analysis finds modes that grow 
at most linearly with time. Choptuik[Q suggests that the presence of the r = 
singularity (persisting even as one takes the M — > and q — > limit) in H results 
in exponentially growing modes rather than the linearly growing modes that we see 
in our analysis, which explicitly assumed regularity. 

8. Constrained Evolution 

In view of Eq. (20), violations of the Hamiltonian constraint lead to spurious h*. 
This can be solved for h*: 

h* = V- 2 R. (47) 

This h* can then be subtracted from the solution at any particular instant, 
reasserting the Hamiltonian constraint. 

Gauge drift, which would arise from an error in setting exactly zero longitudinal 
data, can be similarly suppressed. Section [5] shows how to remove longitudinal data 
from arbitrarily set data. This step can be carried out at any particular instant, 
reasserting the TT requirement. 

9. Summary 

We have shown that gauge and constraint error modes arise in the analysis of the 
3+1 form of the Einstein equations for linearized waves. These modes are shown to 
grow linearly in time and have a form that can be controlled in any well specified 
convergent computational discretization of the evolution equations. By imposing 
the constraints on the free-evolution these modes may be supressed. 
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